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The Cauchy problem for the Liouville equation with a small perturbation 
<9 t 2 w-c^u + 8expw = eF[u]; 0<£«1, (0.1) 

u\ t =o = tpo(x;e), d t u\ t=0 = ipi(x;e), xeR (0.2) 

is considered. We are interesting for asymptotics of the perturbed solution 

u(x, t; e) as e — > 0. 

Perturbation theory for integrable equations remains a very attractive task. 
As a rule a perturbation of smooth solutions such as a single soliton were usually 
considered. We intend here to discuss a perturbation of a singular solution under 
assumption that the perturbed solution has singularities as well. A simple well 
known instance of this kind is a chock wave under weak perturbation as given by 
the Hopf equation with a small perturbation term ut + uu x — ef(u), < e <C 1. 
In this paper we consider a more complicated problem namely a perturbation 
of the singular solution of Liouville equation. We deal with the singularities 
studied by Pogrebkov and Polivanov, [1]. 

1. Unperturbed equation (as e = 0) dfu — d^u + 8expit = was solved 
by Liouville [2] and a formula of the general solution is well known 

r' (s + )r' (s~) 

u(x,t) =ln +V 9 . i~\ \ r = r + (s+)+r-{ S -),s ± =x±t. (1.1) 
r z (s+,s ) 

Initial equations (0.2) give two ODE's for the r± which may be linearized. 

The singularities of the solution occur due to zero of the denominate r(x + 
t,x — t) = r + (x + t) +r~(x — t) = under condition r± > 0. We consider just 
this case. The singular solution is unique under matching condition imposed on 
the singular lines T, [3]. That is continuity (zero jump across the singular lines) 
of two expressions 

[1/2(1% +ul) + 2ex.pu-2u xx ] T = 0, [-u t <t> x + 2u tx ] v = 0. (1.2) 
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2. Perturbed problem (e ^ 0) 

As regards the perturbation operator F[u] we assume that one has no higher 
order derivatives. We desire the singularities do not disperse so that singular 
lines only deform slowly under perturbation. Note that an existence theorem 
for the perturbed problem is not proved up to now. We only give a formal 
asymptotic solution. 

The main goals are both an asymptotic approximation of the singular lines 
and an asymptotic approximation of the solution everywhere out of narrow 
neighborhood of the singular lines. 

Definition 1 The formal asymptotic solution (FAS) of order N is a function 
C/jv(a;,t;e) satisfied to both the equation (0.1) and the initial condition (0.2) up 
to order 0(e N ) everywhere in {x G R, < t < T = const] out of narrow strips 
of order 0(e N ). In the strips there are smooth lines on which the matching 
conditions (1.2) are satisfied. 

Remark. A direct asymptotic expansion as given by u s=a J2 n e ™ ^ i x >t)i £ ' 
does not provide an approach to both the singular lines and the solution inside 
of stripes (near limit singular lines) whose width has the order of O(e). This 
assertion can be verified on a simple example when F = and exact solution is 
taken in the explicit form (1.1) with the smooth functions r±(s ± ;e) depending 
on the parameter e. Asymptotic expansion of such solution as e — > has coeffi- 
cients with the increasing order singularities 0(r(f n ) at the limit singular lines 
r o -"=r + (s+;0) + r_(s-;0) = 0. 

Anzatz. A FAS (for any N) is taken as a finite peace of the asymptotic 
series 

oo 

u(x,t;e) « ^e n u (x,t;e), e -> 0. (2.1) 

The leading order term is here taken as an exact solution of the unperturbed 
equation 

o , r + (s+; e)r'_(s~; e) r . . . 

We permit dependence on a small parameter in the functions r±(s±;e). For 
these functions an asymptotics is constructed in the form 

oo 

r^e^^Wi^), £ ^0. (2.3) 

n=0 

The formulas (2.1)-(2.3) are usually named Bogolubov - Krylov anzatz. This 
approach provides an asymptotic approximation of the singular lines up to any 
order from the equation r + (s ± ; e) + r_(s ± ; e) = 0. Thus the matter is reduced 

to define the coefficients u, r± . 

3. Linearized problem for the correction 
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Corrections u, (n > 1) are determined from the linear equations 
d 2 t u -flg u +8-^ «=/ (z, i; e), {x, t) e i? 2 



with the corresponding initial data. The right sides are here calculated from 

o 1 1 

previous steps, for example /= F[u], /= SF[u] u —Ar' + r'_(u) 2 /r 2 . 

General solution of the homogeneous linear equation is given by the formula 

n 

u(x,t;e) = f + /r' + + j'_/r'_ - 2{j+ +j-)/r. Here j± =j± {s ± ;s) are arbitrary 
functions. In context of the Cauchy problem they are determined by the initial 
functions and an expression similar to the D'Alembert formula takes place. Take 
into account the singularities of order 0(1 /r). 

To solve the nonhomogeneous linear equation the similar formula may be 
used with the functions j±(x,t), depending on two variables, which are defined 
from ODE's in the explicit form. So a solution of the linear equation is deter- 
mined by the integral 

u(x,t) = J J K (s + , s~ , a + , er~ ) / (a + , o~ ) da~ da + 
taken over the characteristic triangle. The kernel K is here expressed by 
K(s+, s -,a+,a-) = — J- r\r + (s+)r^(s-) + r + (a+)r-(a-) 

f ( S +)-r_( S ")] [r+(a + )-r-(a-)]}. 



Let us denote by Uj (x, t; e) a part of the correction corresponding to a solution 
of the nonhomogeneous linear equation. One depends on both the r±(s ± ; e) and 

the j± (s , e), (0 < m < n) and one is determined by means of the quite well 

n n 1 n—l 

specific operator Wi=U [r±,j±, j± ]. It is convenient to identify the singular 
part in these representation 

n n 1 n—l n 1 n—l 

U=V [r±,j±, j± } - (2/r) W [r ± , j±, ]± ]. (3.1) 

Singularities play here a role of secular terms. If we take a FAS in the form 
of the direct expansion, then we find that the singularities became stronger on 
each step. One can guess that this effect is due to the poor approximation 
of the perturbed singular lines. Hence the singularities must to be eliminated 
from the corrections. Just this elimination gives us a good approximation of 
the perturbed singular lines]]. Of course a representation (3.1) with smooth 

n n 

functions V, W is possible just only under some restrictions on the perturbation 

1 One can think these ideas are suitable for another problems with singularities under 
perturbations. 
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operator F[u] . In fact, the more general condition on the F[u] is a representation 

n n 

(3.1) with smooth functions V, W on each step n. 
4. Identification of corrections 

The system of linear equations are based on the leading term whose param- 
eters r± are determined by initial data from equation: 



In 

n=l 



+ l2 e " [3+ l r '+ + 3- A- - ~ { 3+ + 3- ) + "i J = e), (4-1) 

T f r n f n 1 9 / n ' n f \ 

— + E £ " - (i- A-)' - ~ ( i+ - i- ) + 



r" r" 
4- - -F- -2 



n 

Note that there are unknown functions j± apart from r±. 

A naive approach to these equations with a small parameter is as follows. 

Functions r± are identified with the leading term r± =r± (s ± ) which is defined 

n 

from the nonlinear equations as e = 0. After that all functions j± (s ) are 
determined from the recurrent system of linear equations. In this way a direct 
asymptotic expansion is just obtained, whose coefficients have the singularities of 
increasing order 0(r^"") at the limit singular lines ro = r + (s + ; 0)+r_(s~; 0) = 0. 
Deformation of the singular lines is not determined in this way. 

In our approach the functions r± are taken in the form of asymptotic series 
(2.3). Additional ambiguities in the coefficients r± (s ± ), n > 1 are used to 
solve the (4.1), (4. 2) under additional requirements. We desire to eliminate the 

singularities of order C(r _1 ) from the corrections u on each step. Elimination 
of these terms give rise to the algebraic equations 

n n n 1 n— 1 

j+ {s+;s)+j- (s~;s)+ W [r±,j±,..., j±] =0 as r(s+, s"; e) = 0. (4.3) 



Moreover an additional condition is used 



r+ (x)+ r- (x) = 0, Vie R, (n > 1). 



(4.4) 



So we obtain the systems of equations (4.1)-(4.4) for the four functions r±,j± 
on each step n = 1, 2, .... However these equations contain a small parameter 
throw r±(x;e). Hence an asymptotic expansion have to be constructed for the 

n n n,m 

functions j± like (2.3) j± (x;e) « J2 m £m 3± ( x )- Coefficients can be here 
obtained from the recurrent system of algebraic equations 

n.m n,m n,m q 

j+ (s+)+ j- (s~) = W, (n > 1, < m < n) as r (s + ,s") = 0. (4.5) 
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n 71,0 n,m 

So all functions r±,j± (x), (n > 1) and j± (1 < m < n) are denned 
step-by-step from the algebraic and differential equations. The main objects 
are differential equations, which can be represented in the form 

n n 

n' 0' n I > y+ +y- n 

{x)+y-/r- (x) - 2 — ^— - =* (a:), 



71 ' 71' . 71 n , , , , 0,,, 

f n ',° r ^ n '/° V+ -V- l0 (y+ + y-){ r + (x)-r_ (x)) n 
[y + /r + (*)-y_ /r_ (*)) -2^—^+2 ^— =*i (*) 



71 n ™»° 

for the combinations y± (x) =r± (x)+ j ± (x). If we take into account 
r(x, x) —r + (x)+ r°_ (x), then a solution is obtained in the explicit form 

y± = \ j (*o ± .9) r± dx+^ J r± J (v r'/r + (°+ - r_)V r ) dx dx > 

= r(x, a;) J ($fi(x)/r(x,x) + 1 i' a (x)^r + (x) — r_(x)^ /r 2 (x, xj^j dx. 
After that the algebraic equations (4.4)-(4.5) are solved. There are some ar- 

n n.m 

bitrariness in the solution r±, j ± , which however do not effect the obtained 
approximation for the u(x, t; e). 

Conclusion. The main result of the paper is the given above manner of 
determination of the FAS (2.1)-(2.3) in the occurence of singularities. 
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